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Abstract 

It is shown that in the framework of the weak turbulence theory, the autocorrelation and cascade 
timescales are always of the same order of magnitude. This means that, contrary to the general 
belief, any model of turbulence which implies a large number of collisions among wave packets for 
an efficient energy cascade (such as the Iroshnikov-Kraichnan model) are not compatible with the 
weak turbulence theory. 
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I. INTRODUCTION 



Magnetohydrodynamic (MHD) turbulence is present in a wide variety of astrophysical 
systems such as the solar wind, the interstellar medium, accretion discs, and so on. In- 
compressible MHD is the standard model for the study of astrophysical MHD turbulence. 



Although incompressible MHD turbulence has 



decades 
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3een intensively studied for the last several 
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18] many physical aspects 



of the problem still remain unclear. The first model of incompressible MHD turbulence 
was proposed by Iroshnikov [l| and Kraichnan The Iroshnikov-Kraichnan (IK) model 
of MHD turbulence is based on the fact that nonlinear interaction is possible only among 
Alfven waves propagating in opposite directions along the mean magnetic field. Therefore, 
an energy cascade occurs as a result of collisions between oppositely propagating Alfven 
waves. Consider the isotropic excitation of Alfven waves on some outer scale Iq with a char- 
acteristic velocity t> ^ Va, where Va denotes the Alfven velocity. The IK model assumes 
that the energy transfer is local and isotropic in the wave number space. The characteristic 
time scale of the Alfven wave collision is ~ (feV^) -1 , where k is the wave number. Using 
the governing equations of incompressible MHD, it can be shown that during one collision 
the distortion of each wave packet 5vi is of the order Svi/vi ~ vi/Va «C 1. Because these 
perturbations are summed with random phases iV ~ (vi/5vi) 2 ~ {Va/vi) 2 , collisions are 
necessary to achieve the distortion of order unity. Therefore, for the energy cascade time 
r/^f we have 

- 1V * a) 



n-IK 

T ~ 

cas 



kvi vi 

Taking into account the relations e ~ vf/r cas and vf ~ k£k, where e is the energy cascade 
rate and E^ is the one dimensional energy spectrum, we obtain 



£ 



IK 



(eV A ) 1/2 k-^, 



(2) 



which represents the IK spectrum of incompressible MHD turbulence. Due to the condition 
vq Va the IK model is usually deemed as the model of weak MHD turbulence (i.e., 
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20| is 



the model for which perturbation theory called weak turbulence theory (WTT) 
applicable) . 

The IK model for MHD turbulence is isotropic. However, the presence of a mean mag- 
netic field has a strong effect on the turbulence properties, in contrast to a mean flow in 
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hydrodynamic turbulence, which can be eliminated by the corresponding Galilean transfor- 



mation. 



simulations 



?he anisotropy of MHD turbulence had been seen in various numerical 
A 0, Q, 12]. Different theoretical models of weak Q, H, Q as well as strong 2, fl H, Q 
anisotropic MHD turbulence have been developed to the date. In Ref. js| it was shown that 
the lowest order three wave interactions among Alfven waves are possible only if the parallel 
(with respect to the constant magnetic field) wave number of one mode is zero. This implies 
that in the framework of WTT there is no parallel cascade of energy. Consequently, the 
turbulence is anisotropic and the energy is cascaded to larger values of the perpendicular 
wave number. 

Comprehensive study of the MHD turbulence in the framework of the WTT has been 
performed in Ref. {9], where the full coupled equation for shear and pseudo Alfven waves has 
been derived and analyzed. Stationary solutions of the WTT equations has been found and 
was shown that for balanced turbulence the three dimensional energy spectrum Eikn, k±) ~ 
e^-^kj 3 - The validity criterium of the WTT for anisotropic MHD turbulence is js] 

k±vi 

Wa <<: (3) 

The usual interpretation of this condition implies that the nonlinear strain time should be 
less then the wave period. 

In the present paper we study incompressible MHD turbulence in the framework of the 
weak coupling approximation (WCA) [19( , which represents one of the equivalent formula- 
tions of the direct interaction approximation (DIA) |22j. The WCA allows us to study both 
weak and strong limits of MHD turbulence in the framework of a unified formalism. In the 
case of zero residual energy we rederive the WTT equations from the WCA equations of 



incompressible MHD turbulence derived earlier in Ref. [9]. We show that in the framework 
of the WTT the autocorrelation and cascade (evolution) timescales are always of the same 
order of magnitude. This means that, contrary to the general belief, in the WTT N ~ 1 
and, consequently, any model for MHD turbulence that implies N ^> 1 (such as the IK 
model) are incompatible with the WTT. 

The paper is organized as follows. The WCA formalism for incompressible MHD turbu- 
lence is described in Sec. HI] The WTT equations are derived and analyzed in Sec. IIHI The 
conclusions are given in Sec. HVl 
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II. THE WCA EQUATIONS FOR INCOMPRESSIBLE MHD TURBULENCE 



Consider incompressible MHD turbulence in the presence of a constant magnetic field Bo 
directed along z axis. The equations of ideal MHD, governing the evolution of fluctuations 
of the Elsasser variables, are 



d t Vt - V A d z V t = -(Wi ■ V)Ui - Wp (4) 
dtWt + V A d z W 1 = -(Ui • V)Wi - (5) 



where = Ho/y/Anp is the Alfven velocity, Wj = Vj + bj and Ui = Vi — bj are 



fluctuations of the Elsasser variables, vi is the turbulent velocity field, bi = Bi/ y/Anp 
denotes the magnetic field of the fluctuations in velocity units, p is the pressure normalized 
by the density and d t = d/dt. 

Performing a Fourier transform, neglecting pseudo Alfven waves which are known to play 



only a passive role for strongly anisotropic MHD turbulence Jfj], defining the unit polarization 
vector of the shear Alfven waves as = k x z, and introducing the amplitudes of the shear 
Alfven waves as 

W£ = ifa&k, u R = iifaeit, (6) 
Eqs. reduce to the following set of equations 

/oo 
^,201-02^1,2, (7) 
-oo 

/oo 
T 1;2 iP l( f> 2 dTl 2 , ( 8 ) 
-oo 

where Ti j2 = z(ek ■ §k 1 )(k ■ e^) is the matrix element of interaction. 

Ap ply ing; the standard technique of WCA, one can then obtain the following set of equa- 
tions [171 ] 

POO 

C fi + = / |T li2 | 2 r+/ 2 -^ 2 , (9) 



,2 



/oo 
\T lt2 \ 2 I+I 2 dFt 
■oo 

/oo 
|Ti )2 | 2 r+/ 2 -d^ 2 , (io) 
-oo 

and similar equations for and J^. Here k = (k, u), the caret denotes the unit vector, ui 
denotes u ki , u; k = k z VA is the frequency of the Alfven wave, dT\ 2 = c? 4 k 1 (i 4 k 2 5 k; _k 1 -k2' anc ^ 
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<5k-ki-k 2 = _ ki — k 2 ) is the Dirac delta function, 

(</>k0i?) = l£ ^k-k', ("M>&) = ^ k <%-k', (11) 

and 

T t = — 7±- (12) 

Eqs. (JSJ)- (JTUJ) are useless unless some assumptions are made about the frequency de- 
pendence of and 1^. Equivalently, in the framework of the DIA one should make 
some assumptions about the time dependence of G ± (k, r) and Q ± (k, r) [21], which are 
the corresponding inverse Fourier transforms of r^/27r and 1^ with respect to to [i.e., 
(4>k{t + r)d>ti(t)) = Q + (k, r)5 k _ k /] . One of the simplest and frequently used assumptions 



implies 21 



22] 



G ± (k, r) = exp {-\t£\t ± iw)H{r), (13) 

g ± (k, r) = exp {-\Vk\t ± ^ k r)£ k , (14) 

where H(t) is the Heaviside (step) function, and £" k is the energy spectrum. 
Similar to Eqs. (jl3p - (fl4"|) . in the case under consideration we assume 

(£ = Vt (15) 



j± = ^ (in) 

k - (,., -r \2 I („±\2- ^ 



K £ 

Here, for simplicity, we consider the symmetric case rj^ = ry k , and = i? k , which physically 
corresponds to a turbulence with zero cross helicity. Although we consider the symmetric 
case, in the further analysis we will keep the ± signs for the energy spectra in order to 
underline the fact that nonlinear interactions are possible only between counter propagating 
modes. 

r ac = l/?7k is the autocorrelation time. As it was mentioned in Refs. 0, [lj]| the random 
Galilean invariance requires that before applying the WCA (DIA) closure scheme to Eqs. 
(|SJ) one should remove the influence of the velocity field of low frequency modes (for a more 
detailed analysis see Ref. [13]). If this is done, and the corresponding contributions are 
removed from Eqs. fl9|)- f|T0|) . then r ac represents the Lagrangian autocorrelation timescale, 
which is called "the duration of unit act of interaction" in heuristic models of the turbulence. 
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III. DERIVATION AND ANALYSIS OF THE WTT EQUATIONS 



Substituting Eqs. (|T5j) and (fT6j) into Eq. ([9]) and performing an integration with respect 
to the frequencies, we get 



l[UJ 



/oo 
\T lt2 \ 2 R a E+E 2 d!Cl 
-oo 

/oo 
\T h2 \ 2 R b E 2 - <UC k li2 , (17) 
-oo 

/oo 
\T li2 \ 2 R b E 2 -dlCl 2 , (18) 
-oo 



where 



Ra = I ^ + r ? 2 , (19) 

vr fa + r] 2 ) 2 + (u-u} k + Au kl2 ) 2 ' 

^ = 7 \ 77— 7 r , (20) 

{% +VV - 1{ U - U k + Au kl2 ) 

dK\ >2 = d 3 kid 3 k 2 5 k -k 1 -k 2 , and Au kl2 = u; k - u kl + u^ 2 . 

According to Eq. (1161) . 1^ significantly differs from zero when (u — uj^) < i] k . First of all, 
consider the real part of the right hand side of Eq. (fT7|) which describes the nonlinear decay 
of the fluctuations (whereas the real part describes the frequency shift caused by nonlinear 
interactions). Consider the limit 7/k — > 0. In this case, the width of the wave packets tends 
to zero and Eq. ( 1161) yields = E^S^^. In the considered limit the integrals on the right 
hand side of Eq. (TTTT) are dominated by the contribution of the small vicinity of the so-called 
resonant curve (defined by the condition Au kl2 = 0) where the condition 

Au kl2 < r/k, (21) 



is fulfilled. The solution of the resonant condition Auo k \ 2 = is k z2 



3. Consequently, three 



1 0, 3 



The volume 



wave resonant interactions must include the zero frequency mode 
of the wave number space occupied by the resonant area where the condition (I2T1) is ful- 
filled is [r/k/ \du k / <9/c 2 )]£r[. Taking also into account that a typical value of R a and $t(Rb)/iT 
in the resonant area is l/r/k, noting that T 12 ~ k± and assuming that the nonlinear en- 
ergy transfer in the plane is dominated by triad interactions with kj_ ~ fc_|_ 12 then the 
contribution of the resonant area in (say) the first integral of Eq. ( 1171) can be estimated 
as k\E + (k z ,kj_)E~(0,k±)/(duj k /dk z ). Similarly, the contribution of the rest part of the 
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k space is rj\ i k z k\E + {k z , k±)E {k Z) k±)/ujf.. Consequently, the domination of the resonant 
contribution implies 

^E Z ikM 

^ E-(0,k ± ) ^ 1 ) 

If this condition is fulfilled, one can replace R a and ^R,(Rb)/ir by 5(co>k — lj^ + co>k 2 ) = #Acj fel2 - 
Then, the real part of Eq. ( 117]) reduces to 

/oo 
\T ly2 \ 2 E 2 {Ef - E+)5 AuJkl2 d!Ct 2 , (23) 
-00 

where 7k is the total decrement caused by the nonlinear interactions. 

In contrast to the DIA which operates with a two point two time correlation functions 
and/or their Fourier transforms (Q ± (k, r) and J^ 1 ), the WTT implies Markovian closure 
and consequently operates with the Fourier transform of a two point one time correlation 
function n^(t) defined as (0k(£)0k'W) = n k W^k-k'- As known, the nonlinear decrement of 
a one time correlation function for a zero time separation is twice larger then the nonlinear 
decrement of a two point correlation function since all temporal derivatives in the dynamic 
equations now act on both the time variables |2l| . Consequently, according to Eq. (T231) the 
dynamic equation for ri£(t) is 

f) n + r°° 

-gf = ^J |7i, 2 | 2 n 2 -« - <)<W 12 d/C* 2 . (24) 
Similar manipulations lead to the following equation for n^(t) 

~^ = 2nj^ |T lj2 | 2 n+K - n^A^dX*,. (25) 

These equations represent the WTT equations for weak MHD turbulence (with zero residual 
energy), which was first derived in Ref. [9] using the standard WTT technique. 

Let us now turn back to the imaginary part of Eqs. (fT7l) . For validity of the WTT it is 
necessary that the frequency shift caused by nonlinear interactions is smaller than uj^. Taking 
into account that, according to Eq. ( fl8l) . i]^ ~ kj_E~(0, k±)/du)k/dk z , an analysis similar to 
the one performed above shows that, in addition to Eq. ( 1221) . the following condition should 
be satisfied 

— < 1. (26) 
The validity conditions of the WTT in a form similar to Eqs. ([22]) and (|26|) were first 



derived in Ref. 



19]. In general case, the validity criterium of the WTT was found to be 



identical to Eq. (|26|) . The appearance of the additional condition (1221 for the validity of the 
kinetic equation (1211) . is related to the degenerate character of the solution of the resonant 
condition = uo^ — uj^, which implies k^ z — 0, and therefore requires the participation in 
nonlinear interactions of modes with very low frequencies [in the context of the performed 
analysis it is clear that E~(0,kj_) in Eq. (1221) should be understood as the average energy 
density of the modes with k z < r/^/ '(dujk/ 'dk z )}. If Eq. ( 1261) holds but the intensity of the 
low frequency modes is very low, such that left hand side of Eq. (1221 is much greater then 
unity, then the kinetic equation (!!MI) is not valid. However, it can be shown that the WTT is 
still valid then and the nonlinear interactions are dominated by four wave interactions. The 
corresponding kinetic equation was derived in Ref . |J . Introducing the characteristic velocity 
of the perturbations with a characteristic parallel length scale l z ~ l/k z and perpendicular 
length scale l± ~ l/k±, respectively, as vf ~ k z k\E{k Z) k±), the conditions (|22|) and (|26|) 
reduce to the following ones k±vi/k z VA <C 1 and [k±viE(0, k±)]/[k z VAE(k z , k±)] <C 1. Note 
that the first condition coincides with Eq. (J3]). 

Another important temporal characteristic of the turbulence (together with r ac ) is the 
energy cascade timescale T cas which represents the timescale of the energy cascade described 
by Eq. (T24|) . A simple way of determining r cas is the following (a mathematically more sound 
derivation can be found in Ref. 20]): in the case of stationary turbulence, the left hand side 
of Eq. (pH|) is zero. But the characteristic timescale of the energy cascade can be determined 
if we retain only the modes with |kij_| > |kj_| on the right hand side of Eq. (|24l) . The 



equation then obtained describes the energy transfer from the mode with a wave number 
to the modes with higher perpendicular wave numbers. The right hand side can be 
estimated as n + jr cas . The analysis of the left hand side terms similar to the one performed 
for the estimation of 7/k yields 

(27) 

and consequently, the WTT always suggests N ~ ij cas /T ac ) 2 ~ 1. These arguments show 
that the IK model, which implies r ac ~ (/cVa)~\ T cas ~ VA/(kvf) and consequently N ~ 
(Va/vi) 2 3> 1, are incompatible with the WTT. Note that in the considered isotropic case 
Eq. ([3]) leads to the (incorrect) conclusion that the turbulence is weak. On the other hand, 
because T ac kVA ~ 1 Eq. (1261) yields that the IK model does not correspond to the weak 
turbulence limit. 

As was shown above, the width of the resonant are is Ak z ~ (j?k/wk)^ z - Taking into 
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account Eq. (1261) this implies that physically nonlinear interactions in the framework of the 
WTT can be interpreted as resonant interactions among spatially very large wave packets 
[with characteristic size A/ ~ (wk/?|k)/fc z ]. Although the nonlinear interactions are weak, 
the interacting wave packets are very large such that an original wave packet decays before 
interacting wave packets pass through each other. Because the introduction of the Dirac 
delta function in Eq. §M§ requires the limit r/k/^k — ► 0, it is sometimes stated that in the 
framework of the WTT the units of the nonlinear interactions are not wave packets of a 
finite spatial extent but spatially infinite Fourier harmonics 191 ]. 



IV. CONCLUSIONS 



Both strong and weak turbulence regimes of incompressible MHD turbulence were con- 
sidered in the framework of the WCA. We showed that in the framework of the WTT the 
autocorrelation and cascade timescales are always of the same order of magnitude and, con- 
sequently, the WTT always suggests N 3> 1. Physically this is caused by the fact that 
the framework of the WTT can be interpreted as the resonant interaction among spatially 
very large wave packets [with characteristic size Al ~ (ou-k/rfe) /k z ]. Although the nonlinear 
interactions are weak, the interacting wave packets are very large such that an original wave 
packet decays before interacting wave packets pass through each other. 



Acknowledgments 

G.G. and S.M.M. acknowledge the hospitality of the Abdus Salam International Cen- 
ter for Theoretical Physics (ICTP) where part of the work was done. G.G. acknowledges 
partial support from INTAS grant 061000017-9258 and Georgian NSF grants ST06/4-096 
and ST07/4-193. S.M.M.'s work was supported by USDOE Contract No. DE-FG02- 
04ER-54742. These results were obtained in the framework of the projects GOA/2009- 
009 (K.U.Leuven), G.0304.07 (FWO-Vlaanderen) and C 90205 (ESA Prodex 9). Financial 
support by the European Commission through the SOLAIRE Network (MTRN-CT-2006- 



9 



035484) is gratefully acknowledged. 



[1] P. S. Iroshnikov, Astron. Zh. 40, 742 (1963). 
[2] R. H. Kraichnan, Phys. Fluids 8, 1385 (1965). 

[3] J. V. Shebalin, W. H. Matthaeus and D. Montgomery, J. Plasma Phys. 29, 525 (1983). 

[4] J. C. Higdon, Astrophys. J. 285, 109 (1984). 

[5] S. Sridhar and P. Goldreich, Astrophys. J. 432, 612 (1994). 

[6] P. Goldreich and S. Sridhar, Astrophys. J. 438, 763 (1995). 

[7] C. S. Ng, and A. Bhattacharjee, Phys. Plasmas 4, 605 (1997). 

[8] J. Cho and E. Vishniac, Astrophys. J. 539, 273 (2000). 

[9] S. Galtier, S. V. Nazarenko, A. C. Newell, and A. Pouquet, J. Plasma Phys. 63, 447 (2000). 
[10] J. Maron and P. Goldreich, Astrophys. J. 554, 1175 (2001). 
[11] J. Cho, A. Lazarian and E. Vishniac, Astrophys. J. 564, 291 (2002). 
[12] W.-C. Miiller, D. Biskamp and R. Grappin, Phys. Rev. E 67, 066302 (2003). 
[13] W.-C. Miiller and R. Grappin, Phys. Rev. Lett. 95, 114502 (2005). 
[14] S. Oughton and W. H. Matthaeus, Nonlin. Proc. in Geophysics 12, 299 (2005). 
[15] A. Schekochihin, and S. C. Cowley, in Magnetohydrodynamics: Historical Evolution and 

trends, eds. S. Molokov et al., (Springer, 2005). [astro-ph/ 0507686. 
[16] A. Beresnyak and A. Lazarian, Astrophys. J. 640, L175 (2006). 
[17] G. Gogoberidze, Phys. Plasmas 14, 022304 (2007). 
[18] A. Beresnyak and A. Lazarian, Astrophys. J. 682, 1070 (2008). 
[19] B. B. Kadomtsev, Plasma Turbulence (Academic, New York, 1965), p. 54. 
[20] V. E. Zakharov, V. S. L'vov and G. Falkovich, Kolmogorov Spectra of Turbulence I (Springer- 

Verlag, Berlin 1992). 

[21] D. C. Leslie, Developments in the Theory of Turbulence (Oxford University Press, 1973), p. 
24, 93. 

[22] R. H. Kraichnan, J. Fluid Mech. 5, 497 (1959). 
[23] J. A. Krommes, Phys. Reports 360, 1 (2002). 
[24] Y. Kaneda, Phys. Fluids A 5, 2835 (1993). 



10 



